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ABSTRACT 

Aims. We investigate the energy sources of the infrared (IR) emission and their relation to the radio continuum 
emission at various spatial scales within the Scd galaxy M33. 

Methods. We use the data at the Spitzer wavelengths of 24, 70, and 160 /im, as well as recent radio continuum 
maps at 3.6 cm and 20 cm observed with the 100-m Effelsberg telescope and VLA, respectively. We use the wavelet 
transform of these maps to a) separate the diffuse emission components from compact sources, b) compare the 
emission at different wavelengths, and c) study the radio-IR correlation at various spatial scales. An Ha map 
serves as a tracer of the star forming regions and as an indicator of the thermal radio emission. 
Results. The bright HII regions affect the wavelet spectra causing dominant small scales or decreasing trends 
towards the larger scales. The dominant scale of the 70 /im emission is larger than that of the 24 fim emission, 
while the 160 fim emission shows a smooth wavelet spectrum. The radio and Ha maps are well correlated with 
all 3 MIPS maps, although their correlations with the 160 (jm map are weaker. After subtracting the bright HII 
regions, the 24 and 70 /J,m maps show weaker correlations with the 20 cm map than with the 3.6 cm map at most 
scales. We also find a strong correlation between the 3.6 cm and Ha emission at all scales. 

Conclusions. Comparing the results with and without the bright HII regions, we conclude that the IR emission 
is influenced by young, massive stars increasingly with decreasing wavelength from 160 to 24 /im. The radio-IR 
correlations indicate that the warm dust-thermal radio correlation is stronger than the cold dust-nonthermal 
radio correlation at scales smaller than 4kpc. A perfect 3.6cm-Ha correlation implies that extinction has no 
significant effect on Ha emitting structures. 

Key words. Methods: data analysis - ISM: dust - ISM: HII regions - Galaxies: individual: M33 - Infrared: galaxies 
- Radio continuum: galaxies 



1. Introduction 

One of the most important discoveries of the IRAS 
mission was the correlation between the IR and ra- 
dio conti nuum luminosities f o r a sample of galax- 



ies (e.g., Helqu et al. 19851 de Jong et al. 19851: 



Gavazzi et al.lll986f ). iBeck fc Gollal (|l988f) showed 
that this correlation also holds within several spiral 
galaxies including M31 an d M33. The radio-IR cor- 
relation was explain ed by iHelou et all (|l985h and 



de Jong et al.1 (|1985l ) as a direct and linear relation- 
ship between star formation and IR emission. On 
the other hand, there was concern that the cold dust 



emission, as a component of the far-infrared (FIR) 
emission, might not be directly linked to the young 
stellar populat i on. F or the nearby galaxy M31, 



Ho ernes et al.l ([19981 ). using IRAS data, found a 
good correlation between the emission of thermal- 
radio/warm dust and nonthermal-radio/cold dust 
with slightly different slopes. They explained the 
latter correlation by a coupling of the magnetic field 
to the gas mixed with the c old dust. Such a cou- 
pling was also considered bv lNiklas k Beckl (|l997ri 
as the orig i n of the global radio-FIR correlation. 



Hinz et al . (2004) found similar behavior in M33. 
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Table 1. Images of M33 used in this study 



Wavelength Resolution 



20 cm 
20 cm 
3.6 cm 
160 fim 
70 /im 
24 fim 



51" 

540" 

84" 

40" 

18" 

6" 



Telescope 

VLA 1 

Effelsberg 2 

Effelsberg 1 

Spitzer 3 

Spitzer 3 

Spitzer 3 



6570A (Ha) 2" (pixel size) KPN O 4 

1 Tabatabaei et al. (in prep.) 

2 Fletcher. (2001) 

3 Hinz et al. (2004) and this paper 

4 Hoopes fc Walterbos (2000) 



The possibility that the relationship between 
the radio and IR emission might vary within 



galaxies (e.g., iGordon et al. 2004) due to the 



range of individual conditions makes it necessary 
to repeat these kinds of studies with improved 
data in the most nearby galaxies. Furthermore, 
it is uncertain what component of a galaxy 
provides the energy that is absorb ed and re- 
radiated in the IR ( Kennicuttl 1998 ). For exam- 
ple, from the close correspondence between hy- 
drogen re combination line emiss i on and IR mor- 
phologies, Devereux et al. (Il 9941); Devereux et alj 



(}l9960l997l ) and lJones et all (|2002f) argued that the 
FIR is powe r ed predominantly by young O /B stars. 
Deul I dl989l).IWalterbos fc Greenaw"altT (| 19961). and 
Hirashita et al. ( 20031 ) argued that about half of 
the IR emission or more is due to dust heated by a 
diffuse interstellar radiation field that is not dom- 
inat ed by any particula r type of star or star clus- 
ter. ISauvage fe Thuan I (J1992J) suggested that the 
relative role of young stars compared with the dif- 
fuse interstellar radiation field increases with later 
galaxy type. 

Another problem that hinders a better under- 
standing of the radio-IR correlation is the sepa- 
ration of thermal and nonthermal components of 
the radio continuum emission using a constant 
nonthermal spectral index (the standard method). 
Although this assumption leads to a proper estima- 
tion for global studies, it cannot produce a realistic 
image of the nonthermal distribution in highly re- 
solved and local studies within a galaxy, because 
it is unlikely that the nonthermal spectral index 



rema ins constant across a galaxy (|Fletcher et al 
2004) . 

M33 (NGC598) is the neares t late-type spiral 



galax y, at a distance of 840 kpc (jFreedman et al. 
Il99ll) and is ideal to compare the morpholo- 
gies of different components of such a galaxy. 
With an inclination of i = 56° (jRegan fe Vogel 



19941 ). it is seen nearly f ace on. The central posi- 
tion of M33 given by Ide Vaucouleurs fc Leach 
(|l98lh is RA(2000) = l' l 33 m 51.0 s and 
DEC(2000) = 30°39'37.0". The position angle of 
the m ajor axis is PA~ 23° (jDeul fc van der Hulst 
Il987l) . 

With the Multiband Ima ging Photometer 
Spitzer (MIPS. iRieke et al.|[20ol data at 24, 70, 
and 160 ^tm, it is possibl e to in terpret the morphol- 
ogy of M33. iHinz et ah ( 2004 ) compared the first 
epoch data of MIPS with Ha images and radio 
continuum data at 6 cm, using a Fourier filtering 
technique to distinguish different emission compo- 
nents. We now have multiple epochs of MIPS data 
that provide a significantly higher-quality image 
of M33, for example in the suppressing of stripes 
along the direction of scan caused by slow response 
from the far infrared arrays. In addition, w e use a 
wavelet analysis technique that iFrick et al. ( 2001 ) 
have shown to be more robust against noise than 
Fourier filtering. Furthermore, the wavelet tech- 
nique provides more precise and easier analysis of 
the scale distribution of emission energy, especially 
at smaller scales of the emitting structures. We ap- 
ply a 2D-wavelet transformation to separate the 
diffuse emission components from compact sources 
in MIPS mid- and far-infrared (hereafter IR), radio 
(at 3.6 and 20 cm), and Ha images (Table 1). 

The classical correlation between IR and ra- 
dio emission has been at the whole galaxy level. 
However, such a correlation can be misleading 
when a bright, extended central region or an ex- 
tended disk exists in the galactic image. This 
technique gives little information in the case of 



an an ticorrelation on a specific scale (jFrick et al 
2001) . The 'wavelet cross-correlation' introduced 



bv lFrick et alj (|200ll ) calculates the correlation co- 
efficient as a function of the scale of the emitting 
regions. Hughes et al. ( 20061 ) applied this method 



to study the radio-FIR correlation in the Large 
Magellanic Cloud (LMC). 

Hippelein et alJ ( 2003h found evidence for a lo- 



cal radio-FIR correlation in the star forming re- 
gions of M33. In this paper, we investigate this 
correlation not only for the star forming regions 
but also for other structures within M33 at spa- 
tial scales between 0.4 and 4 kpc using the wavelet 
cross-correlation. Instead of the standard method 
to separate the thermal and nonthermal compo- 
nents of the radio continuum emission, we use the 
radio continuum data both at high and low fre- 
quency (3.6 and 20cm, respectively). We compare 
our radio images to the Ha image (which is taken as 
a tracer of the thermal free-free emission) to distin- 
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Fig. 1. The 3.6 cm radio map of M33 (left panel) observed with the 100-m radio telescope in Effelsberg, 
and the combined 20 cm radio map (right panel) from the VLA and Effelsberg observations. The half 
power beam widths (HPBWs) of 84" and 51", respectively, are shown in the lower left corners. 



guish the components of the radio continuum emis- 
sion. 

We study the structural characteristics of the 
MIPS IR images using the 2D-wavelet transfor- 
mation in Sect. 3. The distribution of the emis- 
sion energy at both IR and radio regimes and the 
dominant spatial scale at each wavelength are dis- 
cussed in Sect. 4. We show how the different MIPS 
images are correlated at different scales and how 
the radio-IR correlation varies with components 
of the galaxy (e.g. gas clouds, spiral arms, ex- 
tended central region and extended diffuse emis- 
sion) in Sect. 5. In Sec. 6, we use the Ha map to 
probe the energy sources of IR and radio emission. 
Results and discussion are presented in Sect. 7. An 
ov erview of the pre l imina ry results was presented 



m 



Tabatabaei et all (|2005l ). 



2. Observations and data reduction 

Table 1 summarizes the data used in this work. 
The 3.6 cm radio observations were made with the 
100-m Effelsberg telescop^H during several periods 

1 The 100-m telescope at Effelsberg is operated by 
the Max-Planck-Institut fur Radioastronomie (MPIfR) 
on behalf on the Max-Planck-Gesellscahft. 



(18 nights) between August 2005 and March 2006. 
The reductio n was done in the NOD2 data reduc- 
tion system dHaslamlll974f). The r.m.s. noise after 



combination (|Emerson k. Graevdll988f ) of 36 cov- 
erages is ~220 /uJy/beam. We also observed M33 
at 20 cm with the B-band VLAd D-array during 5 
nights in November 2005 and January 2006 (from 
06 to 08-11-05, 13-11-05, and 06-01-06). The reduc- 
tion, calibration, and mozaicing of 12 fields were 
accomplished using the standard AIPS programs. 
The VLA interferometric data will miss much of 
the extended emission o f the galaxy. For example, 



Viallefond et alj (|1986al ) showed that the WSRT in- 



terferometric map at 1.4 GHz accounted for only 
about 16% of the total emission. Therefore, we com- 
bined tlw_VLA map with the new Effelsberg 20 cm 
map (|FletcherJl200ih to recover the emission from 
extended structures in M33. The r.m.s. noise of the 
final map is ~ 180 /iJy/beam. Detailed descriptions 
of the observations and data reduction of the radio 
data will be given in Tabatabaei et al. (in prep.). 
Both radio maps are shown in Fig. 1. 

2 The VLA is a facility of the National Radio 
Astronomy Observatory. The NRAO is operated by 
Associated Universities, Inc., under contract with the 
National Science Foundation. 
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Fig. 2. Top: the Spitzer MIPS images of M33 at 24 /mi (left panel) and 70 /mi (right panel), 
botto m: the Spitzer MIPS imag e at 160 /mi (left panel) and the Ha map (right panel) from 
KPNO lHoopes fc Walterbod (|200Clh . The resolutions are given in Table 1. 



M33 was ma pped in the IR by MIPS 
( Rieke et al1l2004l ) four times on 29/30 December 
2003, 3 February 2005, 5 September 2005, and 
9/10 January 2006. Each observation consisted of 
medium-rate scan maps with 1/2 array cross-scan 



offsets and covering the full extent of M33. The ba- 
sic data reduction was performed using the MIPS 
instru ment team Data An alysis Tool versions 3.02- 
3.04 (jGordon et al.l 120051 ). At 24/tm extra steps 



were carried out to improve the images includ- 
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Fig. 3. MIPS maps of M33 at 40" resolution (first column) and their wavelet decompositions for 4 different 
scales: 80" (second column), 383" (third column), 626" (forth column), and 1024" (fifth column). At the 
distance of 840 kpc, I" is equivalent to 4pc. The maps at 24, 70, and 160 /xm are shown from top 
to bottom. Before the decomposition, huge HII regions like NGC 604 were subtracted from the original 
images (Sect. 4). Maps are shown in RA-DEC coordinate system and centered at the center of the galaxy. 
The field size is 34' x 41'. 



ing readout offset correction, array averaged back- 
ground subtraction (using a low order polynomial 
fit to each leg, with the region including M33 ex- 
cluded from this fit), and exclusion of the first five 
images in each scan leg due to boost frame tran- 
sients. At 70 and 160 ^m, the extra processing step 
was a pixel dependent background subtraction for 
each map (using a low order polynomial fit, with the 
region including M33 excluded from this fit). The 
background subtraction should not have removed 
real M33 emission as the scan legs are nearly par- 
allel to the minor axis resulting in the background 
regions being far above and below M33. 

The 24 /mi image used consisted of just the 9/10 
January 2006 observations as the depth reached 
in a single mapping was sufficient for this work. 
The 70 /zm image used was a combination of the 
last three observations as the 29/30 December 
2003 70 [im map suffered from significant instru- 
mental residuals. These instrumental residuals were 
much reduced when the operating parameters of 



the 70 /im array were changed after the first M33 
map was made and before the second. The 160 fim 
image used consisted of a combination of all four 
maps. The combination of multiple maps of M33 
taken with different scan mirror angles results in a 
significant suppression of residual instrumental sig- 
natures (seen mainly as low level streaking along 
the scan mirror direction). The images used in this 
work have exposure times of ~100, 120, and 36 sec- 
onds/pixel for 24, 70, and 160 /im, respectively. 

Th e Ha observations by iHoopes fc Walterbod 
(2000) were carried out on the 0.6 meter Burrell- 
Schmidt telescope at the Kitt Peak National 
Observatory, providing a 68' x 68' field of view. The 
MIPS and Ha images are shown in Fig. 2. 

To obtain a proper comparison with the 3.6 cm 
radio map, all images of M33 were convolved to 
the angular resolution of 84". For higher angular 
resolution studies, maps at 18", 40", and 51" were 
also made. The PSF (point spread function) of the 
MIPS data is not Gaussian, in contrast to that of 



6 



Tabatabaei et al.: A multi-scale study of IR and radio emission from M33 



the radio data. Thus, convolutions of the MIPS im- 
ages were made using custom kernels created us- 
ing Fast Fourier Transforms (FFTs) to account for 
the detailed structure of the MIPS PSFs. Details of 
the kernel creation can be found in Gordon et al. 
(2007, in prep.). After convolution, the maps were 
normalized in grid size, rotation and reference co- 
ordinates. Then, they were cut to a common field 
of view (34' x 41' in RA and DEC, respectively). 

In the following sections, we discuss the results 
with and without bright compact sources at each 
resolution to investigate how these sources affect 
the energy distribution at each wavelength and the 
correlations between wavelengths. After convolving 
to each resolution, the bright sources, common to all 
images, were subtracted using Gaussian fits includ- 
ing baselevels ('Ozmapax' program in the NOD2 
data reduction system) . For instance, the giant HII 
complexes NGC604 and NGC595 were subtracted 
at all resolutions. We also subtracted strong steep- 
spectrum background radio sources from the 20 cm 
image. Detailed descriptions of the source subtrac- 
tion at each resolution are given in Sect. 4. 

3. Wavelet analysis of IR emission 

Wavelet analysis is based on a spatial-scale decom- 
position using the convolution of the data with a 
family of self-similar basic functions that depend 
on the scale and location of the structure. Like the 
Fourier transformation, the wavelet transformation 
includes oscillatory functions; however, in the lat- 
ter case these funct ions rapidly decay towards in- 
finity. As a result, iFrick et al.l (|200lh show that 
the wavelet method is more resistant to noise and 
the smoother spectra allow better determination of 
the true frequency structure. The cross-correlation 
of wavelet spectra lets us compare the structures 
of different images systematically as a function of 
scale. 

The wavelet coefficients of a 2D continuous 
wavelet transform are given by: 



1 r +QO x' - x 

W(a,x) = - /(x')V*( )<fa 

a J-co a 



(1) 



where ip(x) is the analysing wavelet, x = (x,y), 
/(x) is a two-dimensional function (an image), and 
a and k are the scale and normalization parame- 
ters, respectively, (the * symbol denotes the com- 
plex conjugate). The above transformation decom- 
poses an image into 'maps' of different scale. In each 
map, structures with the chosen scale are promi- 
nent as they have higher coefficients than those with 
smaller or larger scales. To obtain a good separation 
of scales and to find the scale of dominant structures 
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Fig. 4. The wavelet spectra of the 24 and 70 /xm 
images before (or.) and after (s.s.) subtraction of 
the same sources from the two images at 18" res- 
olution. The data points correspond to the scales 
34, 50, 73, 107, 155, 226, 329, 479, 697, 1015". The 
spectra are shown in arbitrary units. 

in M33, we use the 'Pet-Hat' w avelet that was in- 
troduced bv lFrick et al. (l200lh and applied there 
to NGC 6946. It is defined in Fourier space by the 
formula: 



7T < k < 4tt , 2 - 
7r > k or k > 47r, ^ ' 



where k is the wavevector and k = |k|. 

We decomposed the Spitzer images into 10 dif- 
ferent scales a to compare the morphologies at 24, 
70, and 160 /im at each scal^fl Fig. 3 shows the orig- 
inal maps and the decomposed maps at 4 scales. 
The original 24 and 70 /im maps were smoothed 
to the MIPS 160 PSF with FWHM (full width 
half maximum) = 40" before decomposition. The 
maps at scale a — 80" (0.3 kpc) show the smallest 
detectable emitting structures. Most of the mor- 
phological differences among the MIPS images are 
found at this scale. At scale a = 383" (1.5kpc), the 
prominent structures are spiral arms and the cen- 
ter of the galaxy. The central extended region is 
more pronounced at scale a=628" (2.5 kpc). The 
emission emerges from a diffuse structure at scale 
a =1024" (4 kpc), and it is not possible to distin- 
guish the arm-structure anymore. This structure 
can be identified as an underlying diffuse disk with 
a general radial decrease in intensity. The similarity 
of the 4 kpc maps at different wavelengths indicates 



To have both physically meaningful results and a 
sufficiently large number of independent points, the 
scale a varies between a minimum of about twice the 
resolution and a maximum of about half of the image 
size, a < 1100", for all images studied in this paper. 
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Table 2. Source subtraction thresholds S(A) and number of subtracted sources at different spatial reso- 
lutions. 



Resolution 


Subtracted sources 


S (24 /irn) 


S(70/im) 


S(160/im) 


S(20cm) 


S(3.6cm) 


arcsec 


# 


/ijy/arcsec 2 


/ijy /arcsec 2 


/ijy /arcsec 2 


/ijy/beam 


/ijy/beam 


18 


40 


460 


2650 








51 


15 


115 


1570 


2865 


7190 




84 


11 


85 


900 


2340 


8260 


5780 



that the large scale diffuse emission has the same 
structure at different mid- and far-infrared wave- 
lengths. 

At the smallest scale, the emission emerges from 
spot-like features aligned along filaments with the 
width of the scale. At 24 //m the spots, correspond- 
ing to HII regions, contain most of the energy at this 
scale (see Sect. 6). Diffuse filaments are more pro- 
nounced at 160 /jm . As shown in the next section, 
the fraction of the energy at this scale at 160 is 
less than that at 24 /im. The situation in the 70 (im 
map is in between the 24 and 160 /im maps. It seems 
that the star forming regions provide most of the 
energy of the 24 and 70 /im emission, if the spots 
correspond to these regions, as is discussed in the 
following sections. 

4. Spectral characteristics of IR and radio 
maps 

In this section, we demonstrate how the wavelet 
spectra can be used to investigate the scaling prop- 
erties of the emission. This spectrum is defined as 
the energy in the wavelet coefficients of scale a 
(iFrick et al.ll200lh : 



M(a) 



+ OG /* + 00 



\W(a,x)\ 2 dxL. 



(3) 



After smoothing the 24 /im map to the MIPS 
70 ^m PSF with FWHM = 18", about 40 bright 
sources (mostly corresponding to HII regions) were 
subtracted from both maps (equivalent to remov- 
ing sources with fluxes higher than the lower limits, 
S(A), given in Table 2.). Fig. 4 shows the wavelet 
spectra, M(a), of the 24 /im and 70 /im maps be- 
fore and after source subtraction. Clearly, the small 
scales are the dominant scales before source sub- 
traction. The scale at which the wavelet energy is 
maximum is ~ 70" (280 pc) at 24 /im and ~110" 
(440 pc) at 70 /im. After source subtraction, the 
spectra become flatter. The larger size of the dom- 
inant scale at 70 /im caused by these sources (by a 
factor of ~ 1.6 at this resolution) indicates a slightly 
more extended distribution of dust grains emitting 
at 70 /im than 24 /im in the vicinity of star form- 
ing regions. Moreover, the ratio of the maximum to 



minimum energy in the original spectrum at 24 /im 
is larger (by a factor of 3) than that at 70 /im. This 
indicates that either star forming regions provide 
more energy for the 24 /im emission, or the large- 
scale diffuse emission is stronger at 70 /im than 
24 /im. 

At the next angular resolution, 51", the 
HPBW (half power beam width) of the 20 cm ra- 
dio map, 15 bright sources (HII regions) visible at 
all maps were subtracted from the 20 cm map and 
the smoothed 24, 70, and 160 /im IR maps. In ad- 
dition to these sources, 9 background radio sources 
plus the supernova rem nants SN1, SN2, and SN3 
( Viallefond et al. 1986b) were subtracted from the 
20 cm radio map. The wavelet spectra of the MIPS 
and 20 cm radio maps at this resolution are plotted 
in Fig. 5. The Ha spectrum is also given for com- 
parison. The 24 and 70 /im maps show a smoothed 
version of their distributions in Fig. 4 (the linear 
smoothing factor is ~3 between Figs. 4 and 5). 
However, the effect of the sources can still be seen 
by comparing the left panel with the source sub- 
tracted spectra in the right panel. 

The 160 /im map is hardly influenced by the 
smoothing, as its original resolution is 40" (the 
smoothing factor is ^1.3). The 160 /im spectrum 
is more similar to the 70 /im spectrum than to 
the other spectra. It seems that the compact 
sources have less effect on the energy distribution at 
160 /im, because the smallest scale is not the dom- 
inant scale. Hence, there is no important change in 
the spectrum after the source subtraction. The en- 
ergy shows an increase at the scale of complexes of 
dust and gas clouds (~ 250" or 1 kpc), then a second 
increase in transition to the large-scale structures 
or diffuse dust emission. 

The spectrum of the 20 cm radio image is also 
dominated by bright sources. There is a maximum 
at the scale a = 140", then a decrease towards the 
larger scales with a slope of -0.9. However, a flat 
spectrum remains for the scales less than the size 
of the central extended region (^600") after the 
source subtraction. 
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Fig. 5. The wavelet spectra of the 20 cm radio and IR images at 51" resolution before (left) and after 
(right) subtraction of the same sources. For comparison, the wavelet spectrum of the Ha emission is also 
plotted. The data points correspond to the scales 112, 143, 183, 234, 299, 383, 490, 626, 800, 1024". The 
spectra are shown in arbitrary units. 
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Fig. 6. The wavelet spectra of the 3.6 cm radio and IR images at 84" resolution before (left) and after 
(right) subtraction of the same sources. For comparison, the spectrum of the Ha emission is also plotted. 
The data points correspond to the scales 210, 250, 298, 356, 424, 505, 602, 718, 856, 1020". The spectra 
are shown in arbitrary units. 



The spectra of all the maps at the resolution of 
the 3.6 cm data (HPBW = 84"£3 are shown in Fig. 
6 (left panel). The spectra after subtracting the 11 
brightest HII regions are plotted in the right panel. 
Again, because of these sources, the small scales 
are dominant at 3.6 cm. As shown in comparing the 
left and right panels, they also are seen in all the 
infrared bands, strongly at 24 /im and much more 
weakly at 160 ^m. Here, the smoothing effect is seen 
in all 3 MIPS wavelet spectra. 

Before source subtraction, the 20 and 3.6 cm 
spectra are similar to each other and to the Ha 
spectrum up to the scale of the central extended 

4 Gaussian PSFs with FWHM = 51" and 84" are con- 
sidered to convolve the IR maps to the angular resolu- 
tions of the 20 and 3.6 cm radio maps, respectively. 



region, ^600" (2.5 kpc). All three spectra become 
steeper between a ^600" and a ^850", which 
means that the wavelet energy from structures 
with scales of about half of the size of the galaxy 
(~ 900" or 3.7 kpc) is not as significant as that from 
smaller structures. It seems that a minimum in 
the wavelet spectra of the radio and Ha emission, 
as was also show n in NGC6946 (see Fig. 7 in 
Frick et al.l 12001 ). is a characteristic of this scale. 



However, while this decrease disappears in the 
3.6 cm and Ha spectra after source subtraction, it 
remains in the 20 cm spectrum. This implies that 
besides the bright HII regions there are sources of 
nonthermal emission within the spiral arms and 
central extended region which emitt significantly 
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Table 3. Fractions of the wavelet energy provided 
by the 11 brightest HII regions at different scales 
(Eq. 4). 



A 


A(210") 


A(505") 


A(1020") 


24 fim 


0.94 


0.81 


0.35 


70 /xm 


0.80 


0.57 


0.14 


160 /jm 


0.52 


0.33 


0.11 


3.6 cm 


0.86 


0.69 


0.16 



at 20 cm. 

To estimate how much of the wavelet energy 
M(a) is provided by the subtracted sources, we con- 
sider the following definition: 

M or .(a) - Ms.sXa) 
M or .{a) 



A(a) 



(4) 



where M or .(a) and M s . s .(a) represent the wavelet 
energy before and after source subtraction. Table 
3 shows the fraction of energy produced by the 11 
HII regions in IR and at 3.6 cm for 3 different scales 
(at 84" resolution). The corresponding fractions at 
20 cm are not shown in this table, as some nonther- 
mal radio sources were also subtracted at this wave- 
length. The emission at 24 /mi has the largest en- 
ergy fraction A(a) at all scales. The smallest A(a) 
occurs at 160 /mi. This indicates that HII regions 
play a more important role in providing energy for 
dust emission at 24 /an than at 160 /mi. 

We estimate the uncertainties of the wavelet en- 
ergies of each map by making a noise map with 
the distribution and amplitude of the real noise in 
the corresponding observed map. A linear combi- 
nation of uniform and Gaussian distributions sim- 
ulates the real noise in each observed image. The 
wavelet spectra of the derived noise maps were ob- 
tained using Eq. 3. The resulting wavelet energies 
of the noise maps are taken as the uncertainties of 
the wavelet energies of the observed maps at differ- 
ent scales. The estimated values are at least three 
orders of magnitude less than the wavelet energies 
of the observed maps. Therefore, the results from 
our wavelet spectrum analysis are not substantially 
affected by the noise in the maps. 

5. Wavelet cross-correlations 

A useful method to compare images at different 
wavelengths is the wavelet cross-correlation. The 
wavelet cross-correlation coefficient at scale a is de- 
fined as: 



Aa) = 



J J ITi(a,x) W 2 *(a,x)rfx 
[M 1 (a)M 2 (a)]V2 ' 



(5) 
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Fig. 7. The cross-correlation between 24 and 70 /mi 
(top), 24 and 160 /mi (middle), and 70 and 160 /mi 
(bottom) images at 51" resolution before (or.) and 
after (s.s.) source subtraction. 

where the subscripts refer to two images of the same 
size and linear resolution. The value of r w varies be- 
tween -1 (total anticorrelation) and +1 (total cor- 
relation). Plotting r w against scale shows how well 
structures at different scales are correlated in in- 
tensity and location. The error is estimated by the 
degree of correlation and the number of indepen- 
dent points, n: 



Ar w (a) 



2 

■L\2 



(6) 



where, n = 2.13 (— ) , and L is the size of the maps. 
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First, we examine the cross-correlations be- 
tween the 3 MIPS maps of M33 (Fig.7). There 
are significant correlations between each pair 
of the MIPS maps at diffe rent scales withi n 
0.4<a<4kpc, as r w > 0.75 (jFrick et al.ll200ll ). 



These scales include gas clouds, spiral arms, and 
the central extended region. Comparing the plots, 
the 24-160 /im correlation is weaker than the two 
other correlations, especially at scales smaller than 
the spiral arms (a<400"). This is due to signifi- 
cantly different contributions of HII regions in pro- 
viding energy for dust emission at 24 and 160 /im 
(see Table 3), as the 24-160 /im correlation coeffi- 
cients increase after subtracting these sources. This 
indicates that HII regions influence the correlations 
at scales larger than their sizes (the width of the 
huge HII region NGC604 is about 100" or 0.4 kpc). 

Second, we study the cross-correlations between 
the infrared and the radio continuum emission at 
3.6 cm. As shown in Fig. 8, the correlations be- 
tween 3.6 cm radio and the three infrared emis- 
sion are strong at all scales. At scales smaller than 
the width of the spiral arms (1.6 kpc), the correla- 
tions are higher between 3.6 cm and 24 /im emis- 
sion (and also 70 /mi) than between 3.6 cm and 
160 /im emission. The source subtraction reduces 
the 3.6cm-24/xm and 3.6cm-70/im correlations at 
scales smaller than the spiral arms. This indicates 
that young massive stars are the most common 
and important energy sources of the 3.6 cm, 24 and 
70 /xm emission at these scales. 

Third, the cross-correlations with the radio con- 
tinuum emission at 20 cm are presented in Fig. 9. 
Before the source subtraction, the correlations be- 
tween 20 cm radio and the three IR bands are strong 
at scales smaller than the width of the spiral arms. 
After the source subtraction, the coefficients of the 
IR-20 cm correlations decrease dramatically so that 
they become less than 0.75 at small scales. The rea- 
son is possibly that besides HII regions there are 
other unresolved sources of the 20 cm radio emission 
like supernova remnants that are stronger than at 
3.6 cm but do not emit significantly in the infrared. 
Studies within our Galaxy have demonstrated that 
the ratio of far-infrared to radio continuum emis- 
sion for supernova r emnants is much s maller than 
that for HII regions (jFuerst et al.lll987l ). 

Comparing Fig. 8 with Fig. 9, a drop in cor- 
relation coefficient at the scale of 800" (3.2 kpc) 
is more prominent in the MIPS correlations with 
20 cm than with 3.6 cm. This is due to the strong 
minimum in the 20cm spectrum at this scale (Fig. 
5) that persists even after the source subtraction 
(in contrast to the situation in the 3.6 cm spectrum 
shown in Fig. 6). 



1 .0 r 
0.9 r 
0.8 h 

: 0.7 r 

0.6 h 

0.5 r 
0.4 i_L 

1.1 F 



1 .0 r 
0.9 

□ 0.8 \ 

0.7 I 

0.6 \ 

0.5 \ 
0.4 L 

1.1 fr 

1 .0 r 
0.9 h 
^ 0.8 r 

5 n -i L 

- U . / ; 

0.6 r 

0.5 \ 
0.4 L 



20C 



-<>- 
-A - ~& 



03.6cm-24 / um or. 
A 3.6cm — 24/im s.s. 



-<> <>- 
-A ' 

■-a-' T 



- - A 



O3.6cm-70 / um or. 
A 3.6cm — 70/xm s.s. 



_ 4\' 

- <>--<4> 



A - -A ' - A 



O 3.6cm — 1 60/im or. i 
A 3.6cm — 1 60/im s.S.j 



500 
a (arcsec) 



800 



I 00 



Fig. 8. The cross-correlation between 3.6 cm radio 
and IR before and after subtraction of the same 
sources at 84" resolution. 



6. Comparison with Ha 

If we take the Ha emission as a tracer of star form- 
ing regions that both heat the warm dust and ion- 
ize the gas giving the free-free emission, we expect 
a correlation between Ha emission, IR and thermal 
radio emission. The cross-correlations involving Ha 
and IR images are shown in Fig. 10. The significant 
correlation of Ha with the 24 and 70 /im images at 
the smallest scale confirms that most of the com- 
pact structures in the MIPS decomposed maps at 
24 and 70 /im (Fig. 3) correspond to the HII regions. 
The 160 /im emission is also correlated with Ha 
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Fig. 9. The correlation of the 20 cm radio and IR 
before and after subtraction of the same sources at 
51" resolution (please note that this Fig. starts from 
100" and not 200"). 



emission, although the correlations are weaker than 
the correlations at 24 and 70 fim. It is deduced that 
the role of star forming regions in heating the dust, 
even the cold dust (emitting mainly at 160 /im), is 
generally important at scales smaller than 4kpc. 

Between Ha and 3.6 cm there is a perfect cor- 
relation for both small-scale and extended struc- 
tures: r w >0.95 (bottom panel in Fig. 10). This is 
not surprising as the wavelet spectra of the Ha and 
3.6 cm maps are very similar (see Fig. 6). Thus, not 
only are the star forming regions mostly responsible 
for the 3.6 cm emission, but also the radio contin- 
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Fig. 10. The correlation of Ha with IR (top) at 
51" and with 3.6 and 20 cm radio (bottom) at 84" 
resolution. The maps include all HII regions. 



uum emission is dominated by the thermal free-free 
emission at 3.6 cm for a < 4kpc. 

After smoothing the 20 cm map to the spatial 
resolution of the 3.6 cm map, we obtain the Ha- 
20 cm correlation (Fig. 10). At scales smaller than 
600" or 2.4 kpc, the coefficients are the same as 
for the Ha-3.6cm correlation. This means that the 
bright HII regions are also important sources of the 
20 cm radio emission. It seems that the correlation 
decreases at larger scales, but it is not certain be- 
cause of the large errors. 



7. Discussion 

We studied the spectral characteristics of the MIPS 
mid- and far-infrared, radio and Ha maps and ob- 
tained cross-correlation coefficients between pairs 
of maps. Here, we discuss the implications of our 
results for the understanding of the energy sources 
powering the IR and radio emission from M33. 
Further, we com pare the resu l ts of our analysis on 
M33 with that of lFrick et all <|200lh on NGC 6946. 
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7.1. IR emission 

From the wavelet analysis of the MIPS IR maps the 
following results are found. 

• The 24 and 70 /im emission emerge mostly from 
the bright sources corresponding to star forming 
and HII regions. The influence of HII regions on 
the IR emission is highest at 24 fim (see Table 
3). As expected, at small scales the 24 /im map 
is better correlated with the 70 /im than with 
the 160 /im map. 

• The 160 /im emission emerges from both com- 
pact and extended structures. The compact 
structures correspond to star forming and HII 
regions. As shown in Table 3, the fraction of 
wavelet energy provided by the 11 brightest HII 
regions at 160 /im varies between 1/2 and 1/3 
of that at 24 /zm at different scales. The 160 /im 
map does not show small-scale dominance in its 
spectrum. The 160 /im spectrum shows an in- 
crease when it reaches the scale of complexes of 
dust and gas clouds of ~ 1 kpc, and a second in- 
crease in transition to the large-scale structures 
or diffuse dust emission. 

While it is clear that the ionizing stars heat the 
warm dust, the mechanism heating the cold dust is 
still in debate. According to our wavelet analysis, 
the direct role of the young stellar population de- 
creases with increasing IR wavelength. As shown in 
Table 3, at small scales (< 200" or 0.8 kpc), the rel- 
ative contribution of the 11 brightest HII regions to 
the total emission energy is about 50 % at 160 /im 
which means that the contribution of all HII re- 
gions of M33 is much larger than 50%. Hence, it 
seems that at least up to scales of 0.8 kpc the cold 
dust is effectively heated by UV photons from mas- 
sive ionizing stars. On the other hand, the 160 /im- 
Ha correlation coefficients are smaller than those 
of the 24 /im (or 70/im)-Ha correlation (Fig. 10). 
This means that the number of structures with a 
specific scale emitting strongly in both IR and Ha 
becomes less at 160 /im. This indicates that besides 
UV photons from massive ionizing stars, other heat- 
ing sources contribute also to heating the cold dust. 
These could be non-ionizing UV phot ons either 
from intermediate-mas s stars (5-20 M^) (IXulll99^ 



or from HII regions ( Popescu et al. 2000, 



200 



Misiriotis et al.ll2001 ), or optical photons from sola r 
mass stars ( Xu fc Heloul 1996 : Bianchi et al. 2000h : 
they should contribute to an average radiation field 
as the 160 /im wavelet spectrum is smooth. 

Our findin gs confi r m the argumen ts of e 



(l2002h ;hat the IR is powered predominantly by 
young O/B stars, specifically at 24 and 70 /im. 
This is not necessaril y at variance or in a ccordance 
with the argument of IXu fc HeloiJ (|l996h that only 
27% of the diffuse dust emission from M31 can 
be attributed to the heating by UV photons, be- 
cause of two reasons. First, our wavelet study does 
not include diffuse structures with scales larger 
than 4 kpc. Second, the late-type galaxy M33 is 
more densely populated with young stars than the 
earlier-type galaxy M31. This may mean that the 
relative role of young stars in heating the dust in- 



creas es with later galaxy type (pauvage fc Thuan 



1992J). For a more precise comparison, it is neces- 
sary to accomplish wavelet analysis for other nearby 
gal axies. 



Hinz et al 



Devereux et al 




Devereux et all dl994h : 
1997b and iJones et al 



iHoernes et al.l (|1998l ) showed that the cold dust 
emission in M31 has a weaker correlation with 
the thermal radio than with the nonthermal ra- 
dio emission. Since the contribution of the ther- 
mal emission is higher at 3.6 cm than at 20 cm, 
one may expect a weaker correlation between the 
160 /im and 3.6 cm emission. In M33, we did not find 
significant differences between the 160/im-3.6cm 
and 160/im-20cm correlations (Figs. 8 and 9) at 
scales of 200" (0.8 kpc) < a < 700" (2.8 kpc), even af- 
ter subtracting strong thermal sources (bright HII 
regions) . This could be due to the larger role of UV 
photons from O /B stars in heating the cold dust in 
M33 than in M31, as discussed in the previous para- 
graph. However, a better 160 /im correlation with 
the nonthermal than with the thermal radio is prob- 
able at the scale of the whole galaxy. For instance, 
Hinz et al.l (|2004l ) showed that the morphology of 
the smoothed 160 /im emission from M33 is most 
similar to that of the radio emission at 17.4 cm. 

7.2. Radio emission and radio-IR correlation 

From the wavelet analysis of the radio maps we 
obtain the following results for the 3.6 and 20 cm 
emission from M33. 

• The shape of the 3.6 cm wavelet spectrum is 
very similar to that of the Ha spectrum (Fig. 
6) and, hence, there is a perfect correlation be- 
tween 3.6cm and Ha at all scales (Fig. 10). 
This indicates that not only the star forming re- 
gions are mostly responsible for the radio emis- 
sion at 3.6 cm, but also the thermal (free-free) 
emission is the dominant component at this 
wavelength. But it could also be that the non- 
thermal (synchrotron) emission correlates well 
with the thermal emission. Our recent estima- 
tion of the thermal fraction of the 3.6 cm ra- 
dio continuum emission gives fth ~ 50% for the 
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whole galaxy including the extended disk emis- 
sion (Tabatabaei et al., in prep.). As a higher 
percentage of the thermal fraction is expected at 
scales < 4 kpc, this result is consistent with our 
wavelet study. The nearly parallel 3.6 cm and 
Ha wavelet spectra indicates that extinction is 
scale-independent and has no significant effect 
on Ha emitting structures. 

• The combined 20 cm (VLA+Effelsberg) wavelet 
spectrum represents both compact and diffuse 
structures of thermal and nonthermal emis- 
sion. The bright HII regions (or complexes) give 
strong contribution to the 20 cm emission. The 
strong Ha-20cm correlation at scales smaller 
than the width of the central extended region 
(~ 2.5 kpc) also indicates that the nonthermal 
emission correlates perfectly with the thermal 
emission at these scales. 

• Generally, there are strong 3.6cm-IR correla- 
tions at all scales in the range of 0.8-4 kpc 
within M33 (Fig. 8). The 3.6cm correlations 
with 24 and 70 /iin are better than that with 
160 fim before subtracting bright HII regions. 
However, these correlations become approxi- 
mately the same after subtracting HII regions. 
That is because the 24 and 70 fim emission are 
more influenced by the HII regions than the 
160 /im emission. 

• The strong 20 cm-IR correlations are mainly 
due to the bright HII regions, as the correlations 
fall after source subtraction (Fig. 9). Here, the 
weakest correlation is between the 20 cm and 
the 24 /mi maps. This indicates that the non- 
thermal sources which are associated with star 
forming regions are not as effective as HII re- 
gions in heating the dust. 



In their multi-resolutio n analysis of t he ra dio-IR 
correlation in the LMC. lHughes et alJ (|2006h found 
strong correlations in the regions where the ther- 
mal fraction of the radio emission is high. This 
is in agreement with our results that the 24 and 
70 /im-radio correlations are higher before subtract- 
ing HII regions. Even, the 160/mi-20cm correlation 
follows this pattern before subt racting the therma l 
radio regions. For NGC6946, IFrick et al] (|200l[ ) 
also found a better correlation between the ther- 
mal 3.6 cm radio and IR emission up to scales of the 
width of spiral arms, although th e nonthermal emi s- 
sion dominates in this galaxy (|Ehle fc Becklll993h 



Frick et alJ (|200ll ) for NGC6946, which is also a 
late-type Scd galaxy. 

NGC 6946 was studied at a linear resolution 
of 0.4 kpc corresponding to our analysis at 84" 
HPBW. Comparing Fig. 6 (in this paper) to Fig. 
7 in the Frick et al. ( 200ll ) paper at scales between 
0.8 and 4 kpc, the dominant scale of the Ha wavelet 
spectrum is at the scale of the width of the spi- 
ral arms (1.5 kpc) in NGC 6946, while in M33 it is 
at the smallest scale (^200" or 0.8 kpc in Fig. 6). 
Even after subtracting the 11 brightest HII regions 
the scales corresponding to the spiral arms (~ 400") 
are not dominant in M33. This is because the spi- 
ral arms in M33 are not as pronounced as those in 
NGC 6946. Instead, the Ha emission in M33 has a 
rather clumpy distribution. 

The spectrum of dust emission fro m NGC 6946 
at 12-18 /mi shown in Fig. 7 in the IFrick et al 



( 2001) paper is more similar to the spectrum of dust 
emission from M33 at 160 /im than to that at 24 /im. 
So, the excess of hot dust emission near the massive 
stars found in M33 is not visible in NGC 6946. 

Frick et alJ (|200lh found signatures of 
three-dimensional Kolmogorov-type turbu- 



lence (jKolmogorovl 119411 ) for the ionized gas (Ha 
and thermal radio continuum emission at 3.6 cm) at 
scales less than 0.6 kpc in NGC 6946. An increasing 
part of the spectrum (logl-loga) with a slope of 
5/3 i s typical for this type of turbulence (ppangler 



1991). They obtained this slope considering the 



7.3. Comparison with NGC 6946 

It is instructive to compare the results of the 
wavelet analysis of M33 with those obtained by 



first point at the scale of 0.3 kpc (the resolution of 
their maps was 0.4 kpc), making the slope estimate 
of 5 /3 uncertain. At the same resolution we cannot 
investigate this turbulence for M33, as the smallest 
scale in Fig. 6 is larger than 0.6 kpc. However, at 
higher spatial resolutions the spectra of the radio, 
IR, and Ha emissions from M33 do not show such 
an increase at scales less than 0.6 kpc (Figs. 4 
and 5). The steepest increasing slope is about 1.1 
at scales between 140 and 320 pc (in the 24 /im 
spectrum) (Fig. 4). The slopes become less at 51" 
resolution (Fig. 5). They vary between and 0.6 
at increasing parts of the spectra and between 
and -0.9 at decreasing parts of the spectra. After 
subtracting HII regions, the slopes become even 

flatter. 

Comparing Fig. 12 in the lFrick et al. ( 2001 ) pa- 
per to Fig. 10 (in this paper), the 3.6cm-Ha cross- 
correlation coefficients are less at the small scales 
in NGC 6946 than in M33. As the extinction is pro- 
portional to the dust density and because higher 
dust densities can be found at smaller scales, one 
expects more extinction of Ha at smaller scales. 
This leads to a weaker correlation with the radio 
free-free emission at smaller scales. Hence, it seems 
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that absorption b y dust at scales < 1.2 kpc (Fig. 12, 
Frick et al. 200 ll ) is less significant in M33 than in 
NGC 6946. The existence of considerable extinction 
in NGC 6946 cause d by a pat c hy du st distribution 
has been shown by iTrewhellal (|l998fh 



8. Summary 

Highly resolved and sensitive Spitzer MIPS images 
of M33 at 24, 70, and 160 /^m enabled us to compare 
the morphology of different dust components of this 
galaxy. A 2D-wavelet transformation was used to 
find dominant scales of emitting structures and cor- 
relations of the MIPS images with new radio (3.6 
and 20 cm) maps as well as with an Ha image of 
M33. We compared the results of wavelet analy- 
sis with and without bright HII regions at differ- 
ent resolutions (18", 51", 84"). We found that at a 
characteristic scale of 210" or 0.8 kpc, 80% or more 
of the wavelet energy at 24 and 70 /jm and 3.6 cm 
is from the 11 brightest HII regions. These sources 
cause better correlations between pairs of the 24 
and 70 /im and 3.6 cm maps, while they cause a 
smaller 24-160 /jm correlation. Bright HII regions 
improve the correlations between the 20 cm radio 
map and each of the MIPS IR maps. The Ha emis- 
sion shows a better correlation with IR emission at 
24 and 70 fim than at 160 /im. It is perfectly corre- 
lated with the 3.6 cm radio emission out to scales of 
4 kpc and also with the 20 cm radio emission out to 
scales of the width of the central extended region 
(~ 2.5 kpc). This is understandable because the ra- 
dio continuum emission at 3.6 cm is dominated by 
thermal emission and that at 20 cm by nonthermal 
emission. The most important conclusions of this 
study are: 

• The longer the IR wavelength, the more ex- 
tended is the distribution of dust grains emit- 
ting at that wavelength (Sects. 3 and 4). 

• HII regions influence the IR emission with a 
strength inversely depending on wavelength: 
more influence at 24 /im and less influence at 
160 ^m (Sect. 4). 

• The nonthermal radio emission correlates well 
with the thermal radio and Ha emission out 
to the scale of the central extended region, 
-2.5 kpc (Sect. 7). 

• The warm dust-thermal radio correlation is 
stronger than the cold dust-nonthermal radio 
correlation at scales smaller than 4 kpc (about 
half the size of the galaxy) (Sects. 5 and 7). 

• The effect of extinction in Ha is independent 
of the scale of structures in M33 and is smaller 
than in NGC6946. 
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